Quantum tunneling of a vortex between two pinning potentials 
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A vortex can tunnel between two pinning potentials in an atomic Bose-Einstein condensate on 
a time scale of the order of Is under typical experimental conditions. This makes it possible to 
detect the tunneling experimentally. We calculate the tunneling rate by phenomenologically treating 
vortices as point-like charged particles moving in an inhomogeneous magnetic field. The obtained 
results are in close agreement with numerical simulations based on the stochastic c-field theory. 



Topologically stable quantized vorticity is the smoking 
gun of superfluidity and has been studied extensively in 
superconductors, superfluid Helium and most recently in 
Bose-Einstein condensed (BEC) ultra-cold atomic gases 
[l[ . Tuneability of experimental parameters in dilute-gas 
systems by means of Feshbach resonances and external 
magnetic or electric fields make it possible to investigate 
fundamental condensed matter phenomena in different 
regimes that are not otherwise easily accessible 0, 
Recent experimental advances in the field include demon- 
stration of quantum phase transitions [ij , the Josephson 
effect J5J and vortex formation through synthetic gauge 
fields Q. 

A superfluid vortex in a large BEC is a macroscopic 
object involving many degrees of freedom and is there- 
fore expected to obey the laws of classical physics. The 
static and dynamic vortex properties in a BEC have been 
successfully described by the mean-field Gross-Pitaevskii 
equation [7[ , which is a classical field theory. If the num- 
ber of vortices is so large that their density is compara- 
ble to the particle density, quantum fluctuations domi- 
nate and the mean-field description fails. In this case a 
quantum phase transition is predicted to take place from 
the BEC ground state to a highly correlated state, sim- 
ilar to the Laughlin state in a quantum Hall liquid Q. 
Quantum corrections are also needed if the dynamics of 
a single vortex over distances comparable to its core size 
(the healing length) are considered. 

Quantum tunneling of particles through a barrier, a 
counterintuitive consequence of the superposition princi- 
ple, is known to happen on atomic scales. As pointed out 
by Schrodinger, quantum mechanics admits the superpo- 
sition of macroscopic states as well [9] . An intense search 
for macroscopic systems that show evidence of the effect 
has so far been successful only for strongly-interacting 
condensed- matter systems [13]. Here, we focus on the 
quantum tunneling of a single superfluid vortex between 
two adjacent pinning sites in a weakly-interacting BEC. 
This process involves a small, controllable fraction of the 
total superfluid, thus presenting a route to larger super- 
positions. We employ two methods of study. First, we 
treat the vortex as a charged point particle in an effective 
inhomogeneous magnetic field, exhibited by the pinning 



potentials. We show that a charged particle feels a dou- 
ble well potential in such a field. We then calculate the 
tunneling rate in the double well potential within the 
Heitler-London approximation. Second, we simulate the 
vortex dynamics within the truncated Wigner approxi- 
mation, which takes into account quantum noise on top 
of the mean-field Gross-Pitaevskii equation. The two the- 
ories give closely comparable predictions for the tunnel- 
ing rate, and a much larger rate than found in previous 
studies of vortex tunneling . 

A dilute-gas BEC containing a vortex consists of a su- 
perfluid whirl around a hole in the condensate. In a ho- 
mogeneous system, vortices located at different positions 
are degenerate. In the presence of a localized impurity, 
the degeneracy is lifted. An impurity that repels atoms 
will attract and thus pin the vortex. A vortex can even 
be pinned by a localized impurity in the presence of har- 
monic trapping as studied in detail in Ref. 



12] 



Here we study a vortex in the presence of two equiva- 
lent pinning potentials. First, we employ mean-field the- 
ory to study the possibility of vortex pinning and post- 
pone the consideration of quantum fluctuations. We as- 
sume that the pinning potentials are located at r = ±ro 
and the energy of the condensate is lifted by A[n(r ) + 
n{— i*o)], where A > is the strength of the pinning 
potentials and n(r) is the density of the condensate in 
the absence of the pinning potentials. We approximate 
the vortex located at r = by the mean-field wave 
function ip(r) = ^/nor exp(i9) / y/ 2£ 2 + r 2 Q, such that 
n(r) = |-0(r)| 2 , £ is the healing length and 9 is the polar 
angle of the vector r. If the core of a vortex is located 
at r = xtq, then the extra energy due to the pinning 
potentials is 



AE(x) ~ Xn 



(x - r ) 2 



(x + r ) 



2e + (x-r ) 2 2e + (x + r Q ) 2 



(1) 



We plot x m [ n such that Ai?(x m ; n ) = mm x [AE(x)] in 
Fig. Q] We see that if the separation between pins is 
larger than roughly twice the healing length, then the 
vortex can be pinned by one of the pins. On the other 
hand, if the separation is smaller than that, the vortex is 
located between the two pins. This bifurcation is caused 
by the appearance of an energy barrier between the two 
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FIG. 1. Equilibrium positions x m \ n of a vortex in the presence 
of two pinning potentials separated by the distance do = 2ro. 
Solid line: analytical result based on Eq. ([T]). Open circles: 
results of numerical simulations based on the solutions of Eq. 
(J2j) . Inset: number of particles inside a circle with the radius 
equal to a pinned position of the vortex. Dashed line: result of 
the numerical simulation. Solid line: analytical result based 
of the assumption that the vortex is structureless. 



pins for the vortex to overcome. 

To compare this with more quantitative results, we 
numerically solve the Gross-Pitaevskii equation in 2D Q : 



ft 2 V 2 
2m 



+ V(r)+ 5 |^(r,i)| 2 U(r,i).(2) 



Here V(r) = V± n (r) + V 2 in (r) + V tra . p (r) consists of three 
contributions. The two pinning potentials are of the form 
^pi„ 2 (r) = V(7rA 2 )exp[-(r±r )7A 2 ] with V being 
the strength of a pinning potential, A is its width. The 
trapping potential is Vtrap(r) = mw}r 2 /2, where r is 
the 2D vector. The two-body interaction strength in 2D 
g = V8t: h 2 a/ (ma z ) is expressed via the s-wave scattering 
length a and the axial oscillator length a z = y/h/muj z . 
The axial confinement is assumed much stronger than 
the radial one a± = y/h/muj± ^> a z . We have solved J ~ 

Eq. ([2]) for a stationary state with a pinned vortex. The Here 6(rj — r) = arctan[(yj - y)/(xj - x)] is the phase 



FIG. 2. Effective magnetic field B along the axis through 
the trap center (solid line) and density n (dashed line). The 
two pinning potentials are separated by the distance do = 2£. 
Bo = noh/q, where no is the bulk density at the trap center. 
In this configuration the vortex is pinned classically (cf. Fig. 
[TJ and a double-well potential appears within the effective 
theory. 



0.7/im at half maximum or by heavy atoms in a sp ecies- 
specific optical potential of double-well shape [14fl . The 
vortex can be created by standard techniques |l5| and 
located on one specific pin by adiabatic ramping proce- 
dures. After a given hold time, the vortex location can be 
measured by imaging after a brief expansion phase with- 
out pinning potentials. This would reveal experimentally 
whether tunneling has occurred. 

A vortex in a density gradient feels a Magnus force, 
which is analogous to the Lorentz force on a charged par- 
ticle in a magnetic field. Here, we further develop a phe- 
nomenological approach based on this analogy 3, II | in 
order to estimate the rate of vortex tunneling between 
the pins. We approximate the vortex located at r by the 
Feynman many-body wave-function [l8| 



(ri,...,rjv 



N 

n 



V>o(ri,...,rjv;r). (3) 



healing length is defined as £ = %/ \j2mgn pi where n p is 
the local density of the condensate at a pinning potential 
in the absence of the vortex. We have chosen A = 0.1a_i_, 
g = 0.01ft 2 /m, Vq — 4h 2 /m and the number of particles 
N = f dr\ip(r)\ 2 ss 2 x 10 5 . This corresponds to the 
chemical potential fi — 25huj± and the pinning potential 
amplitude V~o/(irA 2 ) = 5fi. We find that the vortex can 
be pinned by one of the pins if the separation between 
them is larger than a critical value in agreement with the 
previous result. Both are shown in Fig. Q] 

In the recent experiment with 23 Na atoms [l3j ]. 
uj z /2it = 1kHz and ui±/2tt = 20Hz. This gives a± s» 
4.7^m, £ w 1.6/xm and the bulk density n ps 83/im~ 2 w 
5.7n p . Narrow pinning potentials could be realized by 
blue-detuned lasers focused to a beam of diameter of ~ 



of the atom at r, relative to the center of the vortex at 
r. The real function vanishes at the vortex core and 
approaches ^/n far away from it. When the vortex moves, 
the value r changes and the vortex state acquires a phase 
factor. This is similar to the Aharonov-Bohm effect [l9[ 
and motivates the analogy with a charged particle. 

A particle with charge q moving in a magnetic field 
A(r) along a path T acquires an Aharonov-Bohm phase 
A9 = — jr J r A(r)dr, such that its state vector changes 



as — > 



-iA6 



\ip). This can also be calculated as 



A9 = — Im j r dr(tp\V r ip)- When we substitute into this 
expression the Feynman variational state from Eq. (|3|). 
we obtain A9 V = -J d 2 r' J r dr[V r 6>(r' - r)]p(r'; r). Here 
p(r';r) — N J d 2 r 2 ...d 2 r N '4> 2 (r' , r 2 , rjv; r) is the par- 
ticle density at r' of the BEC with a vortex core at r. 



3 



Comparing the expressions A9 and A9 V for the phase 
change, we conclude that we may regard a vortex as a 
charged particle subject to a vector potential 



A(r) = ^J dV[V r 0(r' - r)]p(r'; 



(4) 



We approximate the density of a vortex at r by p(r'\ r) = 
n(r')\r' - r| 2 /[2£ 2 (r') + |r' - r| 2 ] [7], where n(r') is the 
mean-field density without a vortex from Eq. @ and 
£(r') = h/ y/2mgn(r') the local healing length. From 
Eq. §4§ we calculate the magnetic field B = V x A. It is 
directed along the z axis. The magnitude is proportional 
to the convoluted density, blurred on the size scale of the 
vortex core, as see in Fig. [2j We shall study the dynamics 
of the charged particle in this inhomogeneous magnetic 
field. 

A quantum particle with charge q in a constant mag- 
netic field B = Be z and a square box of size L x L has 
energy levels E n = hui c (n + 1/2), where u> c = qB/m fioj ]. 
Each level is highly degenerate. The lowest Landau level 
is characterized by the wave functions 



: exp 



4P 



2/ 2 



(5) 



with I 2 — h/(qB). Degeneracy corresponds to different 
r q , which are separated by |Ar ? | = l 2 /L. Thus, the 
lowest state will be a superposition of many such wave 
functions with different r q . The particle can be found at 
any r q with equal probability. 

It is difficult to deal with this degeneracy numerically. 
A simplified treatment with an effective Hamiltonian can 
be found when the Landau length I is much smaller than 
the healing length £. Formally, we assume delocalized 
Landau wave functions on spatial domains 1Z of the size ~ 
£ 2 where we approximate B(r) as constant, but account 
for the energy dependence of the lowest Landau level on 
larger length scales. Averaging over the domains 1Z then 
yields the effective Hamiltonian 



ft, 



Pi 

2m 



HqB(r) 
2m 



(6) 



Intuitively, this can be understood as a Zeeman Hamil- 
tonian for a particle with orbital magnetic moment p = 
hq/(2m), which is obtained semiclassically for a charged 
particle orbiting at a radius given by the Landau length 
I, with cyclotron frequency uj c . 

The magnetic field depicted in Fig. [5] suggests that Eq. 
© describes a particle moving in a double well potential. 
The tunneling rate of the particle between two wells gives 
the tunneling rate of the vortex between two pins. 

The vortex tunneling rate can be related to a wave- 
function overlap within the Heitler-London approxima- 
tion ll|, |2l|. Let (i — 1,2) be coherent states 
of the vortex pinned by one of the pinning potentials. 
We construct symmetric and antisymmetric states as 



\i> ± ) = (|Vi) ± I V2))/V 2 (! ±1(^1^2)1). The tunneling 
rate is then given by t v sa (^IV^JV^KV'iIV^)!- The 
quantity (-02 l^pij^a) is approximately equal to VqU p . 
The period of the coherent oscillations is T = 2nh/t v 
and thus 



T : 



2ttH 



V n p \{il)i\^2)\' 



(7) 



The states \ipi) are the solutions of Eq. (|6]) with only the 
left and the right pin respectively. The tunneling rate 
shown by circles in Fig. [3] was calculated from Eq. ([7]) 
based on the mean-field density and numerical solutions 
of Eq. ([5]) . In the following we will compare this result 
with simulations in the truncated Wigner approximation. 

Before proceeding, we briefly analyze an approach for 
calculating the tunneling rate of a vortex found in previ- 
ous studies. In Ref. [Hi E3] it was assumed that the 
vortex is a point-like object, i.e. the detailed density 
structure was not taken into account. This amounts to 
considering the effective theory developed in the previ- 
ous section with a constant magnetic field. Indeed, the 
structureless vortex has a healing length which is zero 
and thus p(r' , r) = tiq in Eq. In this case 

A(r) = no- / d 2 v'[V r 6(v' - r)] = ^(e x y - e y x). (8) 

qj 2q 

This corresponds to the constant magnetic field pointing 
in the third direction B = noh/qe z . The coherent state 
of a charged particle moving in this magnetic field is the 
lowest Landau level centered at arbitrary r q and given by 
Eq. ([S]). It is assumed that this degeneracy is broken by 
the pinning potentials and thus r q is the position of one of 
them. Then the overlap of the two states |r g ) and | — r q ) 
is |(— r«j|r 9 )| = exp (— 2nr 2 n )- The distance between two 
pinning potentials should be of the order of I to ensure 
that the overlap has a reasonably large value. On the 
other hand, / is of the order of the average inter-particle 
distance. The latter is much smaller than the healing 
length, i.e. the size of the vortex core, which was ignored 
in the above analysis. Indeed, we have £/Z s» 18. We 
found that the overlap is ~ exp (—No), where N ~ 1CT 3 
for separations above the bifurcation point in Fig. [T] and 
thus the overlap is practically zero. This would give very 
large values for the period of oscillations in Eq. (J7]). In 
contrast, by accounting for the vortex core structure, we 
have found reasonable values for the period. We attribute 
this to two factors: First, the effective number of particles 
involved in the tunneling process No is reduced. This is 
seen from the inset in Fig. Q] and is due to the reduced 
density at and between the pinning potentials, the finite 
vortex core size and a small deviation between the pin 
and vortex positions. Second, the double well felt by 
the vortex becomes shallow for small separation of the 
pins, allowing for much larger overlap of the localized 
states. Thus the exponential behavior found in the above 
analysis breaks down. 
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FIG. 3. The period of quantum oscillations of a vortex be- 
tween two pining potentials as a function of distance between 
them. The open circles represent the result of the truncated 
Wigner approximation. The open squares are the result of 
the effective theory (see text). Inset: One of the trajectories 
of a vortex. Dashed lines indicate the Xmi n ~ ±0.6£ values 
corresponding to do fa 1.8£ in Fig. [1] 



We perform stochastic simulations in the truncated 
Wigner representation to calculate the tunneling rate nu- 
merically. This method enables one to capture many 
quantum features of the system, being increasingly accu- 
rate for short evolution times 



22|. We first find a sta- 



tionary solution of Eq. ((2]) that corresponds to a pinned 
state. This is one of the states for do/£ larger than the 
bifurcation point in Fig. [T] Quantum noise is added 

where 



We use plane waves, 



to this state as V( r ) - ► VK 1 ") 
fj(v) is an orthonormal basis. 
/j(r) = exp(ik J r)/v / V" and ay are complex-gaussian ran- 
dom variables with a*ctj = Sij/2, sampling vacuum fluc- 
tuations in the Wigner representation. M is chosen to 
represent the physical system, while excluding excess vac- 
uum noise [22|. The results should not change signifi- 
cantly upon appreciable change of M. Two different en- 
ergy cut-offs h 2 k 2 M /2m = iOOhu^ and 4QQhuj^ are used 
to determine values for M and define the low-energy c- 
field region. We then propagate this state in real time by 
solving the time-dependent Eq. ((2|). One of the trajecto- 
ries is presented in Fig. El This corresponds to d w 1.8£ 
in Fig. [TJ It is seen that the vortex performs oscillatory 
motion between the two pins with the period T 14/ u±. 
We average the period of such oscillations over 100 tra- 
jectories and present the result in Fig. |3l Our results are 
indistinguishable on the scale of the figure for the two 
cut-offs. We compare the results of the numerical simu- 
lations with the results of the effective theory given by 
Eq. ([7|) in the same figure. The results essentially agree 
in the vicinity of the bifurcation point. We see that the 
period of the vortex oscillations between two pins can be 
of the order of Is if the separation between pins is of the 
order of two-three lengths. 

In conclusion, we have studied the quantum dynam- 



ics of a vortex in the presence of two pinning potentials. 
In particular, we have demonstrated that the vortex may 
quantum mechanically tunnel between the two pins. The 
time scale of the tunneling is achievable in current exper- 
iments if the separation between two pins is of the order 
of a few healing lengths. If realized experimentally, this 
will demonstrate unambiguously the possibility of macro- 
scopic quantum tunneling. 
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